Trap modulation spectroscopy of the Mott-insulator transition in optical lattices 
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We introduce a new technique to probe the properties of an interacting cold atomic gas that 
can be viewed as a dynamical compressibility measurement. We apply this technique to the study 
of the superfluid to Mott insulator quantum phase transition in one and three dimensions for a 
bosonic gas trapped in an optical lattice. Excitations of the system are detected by time-of-flight 
measurements. The experimental data for the one-dimensional case are in good agreement with the 
results of a time-dependent density matrix renormalization group calculation. 

PACS numbers: 03.65.Xp, 03.75.Lm 



It is now possible to study strongly interacting con- 
densed matter systems using cold atoms trapped in op- 
tical lattices 0, Q ■ Recent experiments have led to the 
realization of a Tonks-Girardcau gas in one-dimension 
(ID) 0, 0] and to the detection of the quantum phase 
transition from a superfluid (SF) to a Mott-insulator 
(MI) in bosonic systems [H, [g, 0], to name just two. 
Equilibrium (ground state) properties of the various con- 
densed phases can be extracted from the interference pat- 
tern of the gas after expansion (see for example [H[ ) . Dy- 
namical properties, however, cannot be characterized in 
this way. In order to have access to the excitation spec- 
trum, and therefore to the dynamical properties, addi- 
tional spectroscopic tools such as Bragg spectroscopy Q 
have been designed. Stoferle et al. [6| introduced mod- 
ulation spectroscopy based on a periodic change of the 
optical lattice depth and analyzed the energy absorption 
using time-of-flight techniques. 

In this Letter we introduce a new technique which en- 
ables a direct probing of the phase diagram and excita- 
tion spectrum of an interacting cold atomic gas. Exci- 
tations in the system are produced by periodically mod- 
ulating the strength of an external harmonic potential 
that confines the atoms. The width of the expanded 
atomic cloud in time-of-flight experiments is then used 
as a measure of the increase in energy of the system due 
to the trap modulation. We use this technique to study 
the superfluid to Mott insulator quantum phase transi- 
tion obtained upon increasing the depth of the optical 
lattice. We observe a noticeably different behavior in the 
ID and 3D results. While the 3D results show a sharp 
transition consistent with mean-field predictions, the ID 
results show a broad transition between the superfluid 
and insulating phases. A comparison with exact time- 
dependent DMRG results in ID for a model system al- 
lows us to identify the region of formation of a "wedding- 
cake" structure in the atomic density profile in ID, with 
coexisting superfluid and insulating regions. The latter 



has been recently observed experimentally in 3D . 

The trap modulation spectroscopy presented here can 
be interpreted as a dynamical compressibility measure- 
ment, performed in a regime where the modulation fre- 
quency is much larger than the inter-well tunneling fre- 
quency. This is in contrast to recent proposals of mea- 
suring the compressibility of the system in the sense 
of a volume reduction following an adiabatic compres- 
sion (Toj . as recently performed in a system of interact- 
ing fermions [n]]. We show below that our technique is 
well suited to characterizing dynamical quantities such 
as, e.g., the change of the excitation spectrum across the 
superfluid to Mott insulator quantum phase transition. 
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FIG. 1: (color online) Trap modulation spectroscopy of a Bose 
gas in an optical lattice, a) In the superfluid regime repeated 
compressions of the trap lead to a compression of the Bose gas 
and hence to excitations, whereas in the Mott insulator regime 
(b) the Bose gas becomes incompressible and no excitation 
takes place. In a) and b) the upper diagram represents the 
initial condition before applying the modulation. 

The working principle of our method is illustrated in 
Fig. 1. A Bose-Einstein gas is held in the combined po- 
tential of a harmonic trap of mean frequency to and an op- 
tical lattice of depth Vb . The frequency of the harmonic 
trap is then increased and decreased in a sinusoidal fash- 
ion, with the instantaneous mean trap frequency given 
by 



uj{t) =lo + Slu sin 2 (fit/2), 



(1) 



2 



where 8uj is the pcak-to-peak modulation depth and 
n the modulation frequency. Figure 2 shows results 
of exact time-dependent density-matrix-renormalization- 
group (tDMRG) calculations of this scheme for a model 
system of TV = 15 particles in one dimension. The 
system's dynamics is modeled using the Bose-Hubbard 
(BH) Hamiltonian [12], with a second order Trotter ex- 
pansion of the Hamiltonian and time-steps of 0.01J [HI 
(where J is the inter- well tunneling energy), in which we 
take advantage of the conserved total number of particles 
TV by projecting onto the corresponding subspace. The 
truncated Hilbert space dimension is up to m = 100, 
while the allowed number of particles per site is D = 8. 
The trap modulation Eq.(l) is accounted for by a term 
^muj(i) 2 d 2 L i 2 fii in the BH-Hamiltonian, where m is 
the atomic mass,di the lattice constant and hi the num- 
ber operator on site i 14j. The BH-Hamiltonian de- 



scribes well the microscopic dynamics in the experiment 
(see below), and Fig. 2 shows the salient features of 
the physical picture described above. In the superfiuid 
regime (Fig. 2(a)), after 8 modulation cycles of the trap 
frequency the distribution of atoms inside the lattice is 
markedly different compared to the initial situation as 
the trap modulation has led to a compression of the su- 
perfiuid towards the central lattice sites. In the Mott 
insulator regime, on the other hand, the initial and fi- 
nal distributions are indistinguishable from each other, 
demonstrating that the system has not been excited by 
the repeated trap compressions. 

Our modulation scheme allows us to detect signatures 
of the superfiuid to Mott insulator transition. Figure 3 
(a) shows the results of the ID numerical simulation for 
the relative energy growth rate during the trap modu- 
lation (obtained from a linear fit to the energy increase 
for small t and rescaled by the energy at t = 0). The 
relative growth rate data are plotted as a function of 
the Hubbard parameter U / (2d J) for different occupation 
numbers n at the central lattice site (typically n — 2 in 
our experiments). A weighted average over different oc- 
cupation numbers was also performed in order to take 
into account the contributions of different ID tubes in 
the deep 2D lattice. The figure shows that the relative 
growth rate for TV = 15 and TV = 10 (triangles and open 
squares, respectively) changes rapidly in the parameter 
region 2 < U/(2dJ) < 8. We find that this corresponds 
to the formation of a "wedding-cake" structure in the 
density profile upon increasing U/(2dJ), as in Fig. 2(a) 
and (b). This restructuring of the density profile is ac- 
companied by sharp features in the relative growth rate 
for TV = 15 and TV = 10, which are largely washed out 
by the numerical averaging over the tubes (black dots). 

In order to test these predictions experimentally, we 
created Bose-Einstein condensates of roughly 6 x 10 4 
rubidium-87 atoms inside a crossed dipole trap [3] of 
mean frequency lj ~ 2ix x 80 Hz. A three-dimensional 
optical lattice with lattice constant d^ = 421 nm was 
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FIG. 2: Numerical simulation of the effects of trap modu- 
lation in a one-dimensional lattice with total atom number 
TV = 15 (corresponding to n = 2) for (a) U/(2dJ) = 2 and (b) 
U/(2dJ) = 9. The filled circles represent the atomic distri- 
bution at t = 0, while the open circles show the distribution 
after 8 compression cycles with 8u)/u = 1. 



then adiabatically (within 100 ms) superimposed on the 
condensate. The final depth Vq of the lattice (in the 
range between 6 E Tec and 22 E rcc , where the recoil en- 
ergy E Tec = 2d 2 ) determined the Hubbard parameter 
U/(2dJ), where U is the on-site interaction energy [lj|. 
We realized the superfiuid to Mott insulator transition 
both in ID (d = 1) and 3D (d = 3). For the ID case 
we ramped two of lattice beams to a maximum depth of 
26 E vec (which resulted in a 2D array of one-dimensional 
tubes @, [I3l) and varied the depth Vq of the third lattice, 
whereas for the 3D case all three lattices had the same 
depth. In both cases we performed various tests (vis- 



ibility of the interference pattern [18j, excitation spec- 
trum adiabaticity of the lattice ramps [l^ |) in order 
to identify the critical lattice depth for entering the Mott 
insulating regime and to ensure that the various regimes 
were reached adiabatically, i.e., without exciting the Bose 
gas. 

Once the final lattice depth was reached, the harmonic 
trap frequency was modulated by sinusoidally modulat- 
ing the power of the crossed dipole trap between the ex- 
treme values Pq and P (l + a) at frequency 17 [2(| ■ The 
resulting variation of the trap frequency uj(t) oc y P(t) 
was reasonably sinusoidal for a < 1, with a modulation 
parameter Slo /to = \J\ + a — 1 . After a few milliseconds 
the modulation was switched off and the optical lattices 
were ramped down in 15ms to Vo — AE Tec . The atoms 
were allowed to thermalizc for a further 10 ms, then both 
the lattice and the dipole trap were suddenly switched 
off, and the atoms were imaged after a time-of-flight of 
23.3 ms. For a fixed lattice depth in the superfiuid regime 
we observed a roughly linear increase with modulation 
time in the width of the expanded Bose gas for up to 
10 ms of trap modulation, and this was independent of 
the modulation frequency Q over a wide range (between 
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?s 0.5 kHz and w 4 kHz). These findings are consistent 
with our numerical simulations. 
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FIG. 3: (a) Numerical simulation of trap modulation spec- 
troscopy in ID with 5uj/lo — 0.22 and Q./2n — 1 kHz. Shown 
here are the resufts for boson numbers JV = 15 (open trian- 
gles) and N — 10 (open squares), corresponding to occupa- 
tion numbers n = 2 and n = 1.5 of the central lattice site, 
respectively, as well as a weighted average (filled circles) over 
total atom numbers N = 5, N = 10 and N = 15. The en- 
ergy growth rates are extracted by taking the average of two 
extreme linear fits to the initial growth curve, leading to a 
relative error of around 15 percent, (b) Experimental results 
in ID (open squares) for 6uj/ui = 0.22 and a modulation time 
fmod = 5 ms and 3D (filled squares) with 8uj/lo = 0.41 and 
imod = 10ms. For comparison, the averaged numerical re- 
sults of Fig. 3 (a) are also shown (solid line). The grey bars 
indicate (from left to right) the critical values of U / (2dJ) for 
the formation of the first Mott insulator shell according to 
Quantum Monte Carlo calculations [2lJ and mean-field the- 
ory, and for the second shell according to mean-field theory, 
respectively. 

The experimental results for the ID and 3D cases are 
shown in Fig. 3 (b), where we plot the width a (normal- 
ized to the width o~q in the absence of any modulation) of 
the Bose gas in time-of-fiight (which reflects the increase 
in energy of the system) after a fixed modulation time 
imod as a function of the Hubbard parameter U/(2dJ) 
for d = 1 and d = 3. In the one-dimensional case, the ex- 
perimental results are consistent with the (averaged) ID 
numerical simulations shown in Fig. 3 (a) if an overall 
scaling factor is applied to the theoretical results. Gener- 
ally, we see a decrease in the excitability of the Bose gas 
reflected by <t/<t approaching unity. While in ID this 
decrease is gradual and even for large values of U/(2dJ) 
there is a residual excitation of the Bose gas, in 3D we 
observe a sharp drop in ct/cto around U/(2dJ) = 4. For 
U/(2dJ) > 5, cr/co remains constant around 1, indicat- 
ing that the system is incompressible. This agrees well 
with the physical picture of a transition from a compress- 
ible (and hence excitable) superfluid to an incompressible 
Mott insulator for U/(2dJ) ~ 5.8 (mean-field prediction) 
or U / (2d J) ?s 4.9 (Quantum Monte-Carlo simulation) 0. 



Performing the above experiment in 3D for different 
values of the modulation depth 5uj/lu, we found that 
the sharp feature around U / (2dJ) = 4 was clearly vis- 
ible for 5lo/lu < 0.5 but became increasingly washed 
out for larger values. Indeed, for an intermediate value 
Suj/lu = 0.7 (see Fig. 4 (a)) we saw a structure with two 
apparent "steps" at U/(2dJ) « 4 and U/(2dJ) « 9, and 
ct/cto approached unity for U / {2d J) > 10. The proximity 
of the step positions to the theoretical values of U/ (2d J) 
for the formation of Mott insulator shells with n = 1 and 
n = 2 atoms per lattice site, respectively, 0, [22j suggests 
the following intuitive explanation. For small modulation 
depths, the formation of the n = 1-shcll greatly reduces 
the overall excitability of the system, while for larger val- 
ues of Suj/lu the n = 1 shell, being located away from the 
trap center and hence exposed to larger potential gradi- 
ents when the trap is compressed, remains partly com- 
pressible and only when the n = 2 shell is formed does 
the system become insensitive to the trap modulation. 

We tested this hypothesis by measuring the Bose gas 
response as a function of Slu/uj in the superfluid and in 
the Mott insulator regime. Figure 4 (b) shows that while 
in the superfluid regime ct/cto increased roughly linearly 
with the modulation depth, in the Mott insulator regime 
we observed that for Suj/lu < 0.6 there was no dependence 
on the modulation depth. For larger values of Suj/uj, the 
Bose gas could again be excited, suggesting that in this 
regime the trap compression was sufficient to overcome 
the potential barrier U for moving atoms between adja- 
cent lattice sites. Alternatively, a strong trap modula- 
tion might lead to a heating of the thermal fraction of 
the sample, thus "melting" the Mott insulator and 
increasing the compressibility of the Bose gas . While we 
were unable to rule out cither of these two scenarios, this 
remains an interesting question for future experiments 
and theoretical analysis. 

Finally, we investigated the response of our system to 
trap modulation for adiabatic and non-adiabatic loading 
of the lattice. In all of the experiments described above, 
we made sure that the timescales and shapes of the in- 
tensity ramps for the lattice beams were such that we 
reached the final lattice depth adiabatically and hence 
always prepared the system in its ground state. One 
can, however, prepare the same final lattice depth non- 
adiabatically, in which case one expects to find the sys- 
tem in an excited state. If the lattice depth is above the 
critical value for the Mott insulator transition, this means 
that for non-adiabatic loading the system will not be in 
a pure Mott insulator state and should, therefore, retain 
its compressibility. We tested this assumption by using 
a two-part ramp for the optical lattice with a 100 ms ex- 
ponential ramp up to 5 E lec and a linear ramp in 0.1 ms 
to the final value Vq. The second part of the ramp was 
highly non-adiabatic, which we verified by lowering the 
lattice after a holding time of 5 ms and checking that the 
Bose gas was excited, as expected. We then performed 
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we obtained good agreement with exact numerical calcu- 
lations. For the three-dimensional case a full theoretical 
interpretation will require further calculations. As a per- 
spective for future studies we believe that it is of great 
importance to explore in detail the differences and sim- 
ilarities of our approach with other methods based on a 
quasi-static compression of the trap 11 1 . 

Financial support by MIUR PRIN-2007 and the E.U. 
projects SCALA and NAMEQUAM is greatfully ac- 
knowledged. Numerical simulations in this work were 
done using the DMRG code released within the PwP 
project (www.dmrg.it). 



FIG. 4: (a) Trap modulation spectroscopy in 3D for Suj/lu — 
0.22 (filled squares) and Su/u = 0.7 (open triangles). The 
grey bars are defined in Fig. 3. (b) The response of the Bose 
gas to trap modulation as a function of the modulation depth 
Suj/lu for U/{2dJ) = 2 (filled squares) and U/(2dJ) = 9 (open 
circles) . 



trap modulation spectroscopy using a short modulation 
time of 5 ms in order to minimize the role of the spuri- 
ous effects on a/ao due to the loading process. Figure 5 
shows that for a non-adiabatic loading ramp, the Bose 
gas does, indeed, remain excitable well above the critical 
value for the formation of a Mott insulator. 
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FIG. 5: Trap modulation spectroscopy for an adiabatically 
loaded lattice (filled squares) and for non-adiabatic loading 
(open circles). In this experiment, Soj/uj = 0.41, and the grey 
bars are defined in Fig. 3. 



In summary, we used trap modulation spectroscopy 
in order to study the superfluid to Mott-insulator phase 
transition. The method presented here seems to be a 
promising tool for the characterization of strongly inter- 
acting cold atomic systems. For the one-dimensional case 
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